Relativistic corrections and non-Gaussianity in radio continuum surveys 
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Forthcoming radio continuum surveys will cover large volumes of the observable Universe and will 
reach to high redshifts, making them potentially powerful probes of dark energy, modified gravity 
and non-Gaussianity. We consider the continuum surveys with LOFAR, WSRT and ASKAP, and 
examples of continuum surveys with the SKA. We extend recent work on these surveys by including 
redshift space distortions and lensing convergence in the radio source auto-correlation. In addition 
we compute the general relativistic (GR) corrections to the angular power spectrum. These GR 
corrections to the standard Newtonian analysis of the power spectrum become significant on scales 
near and beyond the Hubble scale at each redshift. We find that the GR corrections are at most 
percent-level in LOFAR, WODAN and EMU surveys, but they can produce 0(10%) changes for 
high enough sensitivity SKA continuum surveys. The signal is however dominated by cosmic vari- 
ance, and multiple-tracer techniques will be needed to overcome this problem. The GR corrections 
are suppressed in continuum surveys because of the integration over redshift - we expect that GR 
corrections will be enhanced for future SKA HI surveys in which the source redshifts will be known. 
We also provide predictions for the angular power spectra in the case where the primordial pertur- 
bations have local non-Gaussianity. We find that non-Gaussianity dominates over GR corrections, 
and rises above cosmic variance when /nl > 5 for SKA continuum surveys. 



PACS numbers: 



I. INTRODUCTION 



Radio continuum surveys for cosmology are entering 
a new phase, given the imminent surveys with LOFAR 
(the LOw Frequency ARray for radio astronomy [1]), 
WSRT (Westerbork Synthesis Radio Telescope [2]) and 
ASKAP (Australian SKA Pathfinder [3]) telescopes, and 
the prospect of the Square Kilometre Array (SKA) in the 
coming decade. Increased sensitivity, a very wide sky cov- 
erage, and deep redshift reach will facilitate cosmological 
observations with significant accuracy. 

This has only recently been explored in [4], which 
analyzed what can be achieved by surveys with LO- 
FAR [5], WSRT (WODAN, Westerbork Observations of 
the Deep Apcrtif Northern sky survey [6]) and ASKAP 
(EMU, Evolutionary Map of the Universe [7]), via three 
experiments: auto-correlation of radio sources, cross- 
correlation of radio sources with the Cosmic Microwave 
Background (the late Integrated Sachs- Wolfe effect), and 
cross-correlation of radio sources with foreground ob- 
jects (cosmic magnification). The auto-correlation func- 
tion has been further investigated by [8] , which examines 
the impact of cross-identification of radio sources with 
optical redshift surveys. 

The huge volumes covered by forthcoming radio sur- 
veys, and their deep redshift reach in comparison to 
current and future optical surveys, mean that correla- 
tions on scales above the Hubble horizon H~ 1 (z) will 
be measured. On these scales, the standard analysis of 
the power spectrum is inadequate - because this anal- 



ysis is Newtonian, i.e. it is based on the assumption of 
sub-Hubble scales. The Newtonian analysis must be re- 
placed by the correct general relativistic (GR) analysis 
in order to consistently incorporate super-Hubble scales. 
On small scales, the Newtonian analysis is a very good 
approximation, but on larger and larger scales, the GR 
corrections become more significant. 

On these larger scales, any primordial non-Gaussianity 
in the matter distribution also grows larger. This probe of 
non-Gaussianity is expected to become competitive with 
the CMB for large-volume surveys such as those in the 
radio. Thus it is important to perform a GR analysis in 
order to correctly identify the non-Gaussian signal. 

Unfortunately cosmic variance also becomes more and 
more of a problem on these larger scales covered by radio 
surveys. However, it is possible to beat down cosmic vari- 
ance by using different tracers of the underlying matter 
distribution. 



In this paper, we re- visit the Newtonian analysis of ra- 
dio continuum surveys, and include for the first time the 
terms that were not considered in the auto-correlations 
computed by [4, 8] - i.e. redshift space distortions, lens- 
ing convergence and the GR corrections (potential and 
velocity terms). We do this first in the Gaussian case, 
and then when there is primordial non-Gaussianity of 
the local type. 
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II. GENERAL RELATIVISTIC CORRECTIONS 
TO THE ANGULAR POWER SPECTRUM 

A GR analysis of the matter power spectrum needs to 
start by correctly identifying the galaxy overdensity A 
that is observed on the past light cone. In the standard 
Newtonian approach, the overdensity S is defined in some 
gauge. A change of gauge gives effectively the same re- 
sults on sub-Hubble scales, but leads to different results 
on large scales - and this remains true even if we use 
gauge-invariant definitions of 5. There are many gauge- 
invariant definitions of the overdensity, but the observed 
A is unique, and is necessarily gauge- invariant. In addi- 
tion, we need to account for the distortions arising from 
observing on the past light cone, including all redshift 
space effects and volume distortions. 

The GR analysis of the power spectrum has recently 
been developed [9-21], and the consequences for the cor- 
relation function at wide angles and large scales have also 
been computed [22]. The observed overdensity in direc- 
tion n at redshift z is 



A(n,z) = 8 z {n,z) 



5V(n,z) 



(1) 



where S z is the redshift space density perturbation and 
V is the physical volume density per redshift interval 
per solid angle [11]. Each term on the right is physically 
defined and hence gauge-invariant. In order to compute 
these terms, we can choose any gauge that we prefer. 
Here we work in the Newtonian gauge, 



ds 2 = a(rj) 2 - (1 + 2$)dr) 2 + (1 - 2^)dx' 



(2) 



In the Newtonian approach, the standard observed 
overdensity includes redshift space distortion and lens- 
ing contributions [23]. Lensing is often omitted in the 
Newtonian approach, although the lensing contribution 
to the density contrast can be significant [24]. In their 
analysis of radio continuum surveys, [4] did not include 
either the redshift space distortions or lensing - and our 
angular power spectra here generalize their results also 
within the standard Newtonian approximation. The stan- 
dard observed overdensity is 



1 dv 



A = S - -n ■ — + (5s - 2)k, 
H dx 



(3) 



where \ is the comoving radial position of the source, v 
is the peculiar velocity, and n is the lensing convergence, 

- = ~K f dfj ^ ($ + n (4) 

The slope s characterizes the change of the number den- 
sity with respect to the threshold magnitude [12]: 

d 2 
s = -^-\og 10 N(z,m < to*) = - ((a - 1) + 1), (5) 



where the parameter (a — 1) is used in [4]. 

The observed overdensity in GR is given by [11, 12] 

^GR. ^std _|_ ^GR corr ^g^ 

A GRcorr = + 1)$ + ( 5s _ 2 )$ + 1$' 

+ A / dr](& + + i -i- / d77($ + *) 



An ■ v, 



where 



, W 2 -5s r 



(7) 



(8) 



The GR correction includes potential contributions (both 
local at the source and integrated along the line of sight) 
and a Doppler velocity contribution. 

A further gauge subtlety arises in relation to the bias 
of the sources. The simple bias relation <5 = b5 m , where b 
depends only on redshift and not on scale, may be applied 
in any gauge on sub-Hubble scales. But on larger scales, 
the relation is gauge-dependent, and we need a GR analy- 
sis. The simple bias law is applicable in the synchronous- 
comoving gauge, as shown by analysis of the spherical 
collapse model and by the physical argument that galax- 
ies and dark matter follow the same velocity field [11-14]. 
Then we must transform from synchronous-comoving in 
order to get the correct bias relation in Newtonian gauge: 



III. 



aH ,, 
3— (6-l>. 



OBSERVABLE EFFECTS OF GR 
CORRECTIONS 



(9) 



In a radio continuum survey, sources are detected via a 
standard approach, as used in [4, 8]. First the imaging is 
done (using e.g. CLEAN, maximum entropy, or sparsity 
methods), and then discrete objects are detected using 
SExtractor or similar codes. (This is to be distinguished 
from intensity mapping.) The sources that are isolated 
from images generated from the survey are used to com- 
pute the power spectrum. Since the surveys integrate over 
redshift, the overdensity in a given direction may be ex- 
panded in spherical harmonics as 



A(n) = y^ y at m Yt m (n), 



(10) 



(if,, 



dn n dzY; m {n)W(z)A(n,z), (11) 



where W(z) is the number count per unit solid angle at 
redshift z. Then the integrated angular power spectrum 
is 



C, 



(|a& 



(12) 



= - d\nk dz dz' V in (k)F e (k, z)W(z)F£ (k, z')W(z'), 
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FIG. 1: Fractional error ACf /Cf R defined in (14), when the various types of terms X in the full result Cf n are neglected. 
This is shown for the LOFAR MS 3 , EMU and WODAN continuum surveys. Here and in Fig. 4, a broken curve denotes the 
absolute value of a negative quantity. The GR correction for LOFAR is shown on its own in Fig. 4. 
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FIG. 2: Number density for SKA-like continuum survey with 
different sensitivities S, in units (arcmin)" 2 . (The number 
counts W in (11) are per solid angle.) 



Here the f a incorporate the auto- and cross-correlations 
of the various terms in (3) and (7). 

We can compute (12) with and without GR correc- 
tions in order to assess their impact. It is useful to define 
(following [11, 12]) 



AC} 



r<GR 



C 



GR-X 



(14) 



where GR denotes the Ci calculated with the full GR 
overdensity (6) and X denotes the contribution from var- 
ious terms in (3) and (7): 

X = redshift term (ex dn ■ v /dx) , 
tensing term (oc k), 
velocity term (oc n ■ v), 
potential terms (all those involving <£>, 
GR correction term A GRcorr . 

Therefore ACf gives the contribution of X in each case. 
A negative AC^ corr means that GR corrections reduce 
C e . 



where Tin is the (dimensionless) initial power spectrum LOFAR, EMU, WODAN 

of $, and the function Fg is of the form [11] 

We adopt the models of [4] for radio source number 
F e (k,z) = je(kx)fi(k,z) + j' e {kx)f2(k,z) density n(z) and bias b(z) for the LOFAR, EMU and 

WODAN surveys and we use their luminosity function 
parameter (a — 1)(= 2.5(s — 1)) for each survey. (See 



+ f/(kx)f3(k,z)+ / dxj t (kx)U(k,z). (13) 
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FIG. 3: The same as Fig. 1, for future SKA-like continuum surveys with different flux limits S. For the highest sensitivity, the 
GR correction is shown on its own in Fig. 4. 



Fig. 1 in [4] for n(z).) For illustrative purposes, we use 
the higher flux density threshold LOFAR MS 3 number 
densities; the LOFAR Tier 1 results are similar to the 
EMU results. Assuming Gaussian perturbations, for the 
surveys studied in [4], the various terms in (3) and (7) 
are shown in Fig. 1. The curves are the absolute value 
of AC^ defined in (14), rescaled by the spectra with full 
GR corrections C^* R . 

The GR corrections are percent-level at low multipoles 
I < few, and are largest in the LOFAR MS 3 case. 

SKA 

We obtain models of n(z) for SKA-like continuum 
surveys (i.e. without redshifts) from the S 3 simulations 
of [25]. We examine both moderate and ultra-deep SKA- 
like surveys, with flux density thresholds of 1/J.Jy, 0.1 
/iJy and O.Ol^Jy, applying the bias models of [4] for 
each source type. This highest sensitivity survey will be 
costly to achieve, but we examine it to see whether GR 
corrections are important in the extreme case. In ad- 
dition, we measure the luminosity function parameter 
(a — 1)(= 2.5(s — 1)) using the S 3 simulations. We find 
(q — 1) of —0.20, —0.40 and —0.55 (corresponding to a 
slope s = 0.32,0.24 and 0.18) for the 1/xJy, O.LuJy and 
0.01/iJy surveys respectively. The number density for var- 
ious sensitivities is shown in Fig. 2. 

The resulting size of the various contributions and of 
the GR corrections is shown in Fig. 3. For the highest 



sensitivity SKA, the GR effects are at the O(10)% level 
for £ < few. 



Cosmic variance 

Figure 4 shows the GR correction signal within the cos- 
mic variance (shaded bands): for the LOFAR MS 3 (solid 
black in the left panel) and for the highest sensitivity 
ultra-SKA (S > 0.01/iJy) case where the GR correction 
is most visible (solid black in the right panel). The cos- 
mic variance for a survey covering a fraction / s k y of the 
sky is 

° {C ' ] ~ ^rW°- < 15 » 

This is highest on the largest scales, where the GR cor- 
rections are strongest. It is not surprising that cosmic 
variance overwhelms the GR correction signal, as shown 
in Fig. 4. However, cosmic variance can be overcome by 
using multiple tracers of the underlying matter distri- 
bution [26-28]. In principle the GR corrections in radio 
continuum surveys may be detectable using this method. 
Furthermore, if the observations include source redshifts 
- i.e. if we have an HI survey - then we expect that the 
GR corrections will be stronger than in the continuum 
case, where redshift is integrated over. 
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FIG. 4: The fractional GR corrections, IAQ 



r (/NL)/Cf R (/ NL )|, for LOFAR MS 3 and ultra-SKA (S > O.Ol^tJy) continuum 



surveys, including cases with non-Gaussianity. Cosmic variance is shown by shaded regions. (Note that the correction is negative 
in all cases.) 



IV. NON-GAUSSIANITY AND GR 
CORRECTIONS 

Several theories of the very early universe predict non- 
Gaussian initial conditions for the probability distribu- 
tion function of cosmological perturbations. One way to 
test this is via the large-scale clustering of matter. Among 
the numerous models of non-Gaussianity considered (sec 
[34] and references therein), one of the most studied is 
the so-called local form, in which deviations from Gaus- 
sianity are parametrized as: 

$ = s -/nl(^-<^>), ( 16 ) 

where <f> g is Gaussian. This leads to a scale-dependent 
correction to the bias [13, 29-33]: 

b(k, Z ) = b g {z) + y NL [b g (z) - 1] J^)d\ z) . ( 17 ) 

where b g is the Gaussian bias, <5 ec is the critical matter 
ovcrdcnsity for ellipsoidal collapse, T{k) is the transfer 
function, D{z) is the growing mode of density perturba- 
tions, and c is the speed of light. Following [32], we take 
<J ec = 1.68^0^75= 1-45. 



In a radio continuum survey we count the number of 
galaxies (as opposed to intensity mapping), and so we 
expect that galaxy number density peaks are related to 
the high density peaks in the overall matter distribution 
in the usual way. Therefore it should be reasonable to 
apply (17). 

A scale- independent bias with / NL = 0, as in (9), be- 
comes scale-dependent when /nl 0, by (17). Therefore 
the non-Gaussianity can in principle change the shape 
of the angular power spectrum on large scales. The non- 
Gaussian signal clearly grows with scale - as does the 
GR correction. In order to obtain the correct prediction 
of non-Gaussianity, it is therefore necessary to use the 
full GR power spectrum, as pointed out in [13]. Although 
the GR corrections arc independent of non-Gaussianity, 
the full Cf K includes cross-correlations of GR correction 
terms with non-Gaussian correction terms, so that the 
difference between Cf R and C| td can vary with /nl- 

Figure 5 shows the predicted Ct with various values 
of /nl for different radio surveys. We have included full 
GR corrections in all spectra. We can see a clear deviation 
on large scales that overcomes cosmic variance for a non- 
Gaussianity signal of /nl ^ 20 in the LOFAR MS 3 case, 
and /nl ^ 5 for a SKA continuum survey at 1/xJy or 
fainter. 
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FIG. 5: The imprint of primordial non-Gaussianity on the observed galaxy-galaxy angular auto-correlation spectra for various 
radio continuum surveys, including full GR corrections. Shaded bands show the cosmic variance. 



This can also be seen in Fig. 6, where the fractional 
differences between the cases with and without non- 
Gaussianity are plotted. Comparing the spectra with 
(solid) and without (dashed) GR corrections, it is clear 
that ignoring the GR correction has a marginal effect on 
the detectable non-Gaussian signal for the sensitivity of 
LOFAR MS 3 , and a small effect (though growing with 
increasing /nl and scale) on ultra- SKA (S > 0.01/iJy). 

The presence of dominant non-Gaussianity has a size- 
able effect on the detectability of GR corrections. As 
Fig. 4 shows, in most cases a large /nl makes the GR 
corrections relatively less important. This is understand- 
able since /nl only boosts the 'standard' term (3) (via 
the overdensity term 5), which makes the GR correction 
terms less visible. 

It is however worth pointing out that in radio surveys 
the radial galaxy distribution and the bias are particu- 
larly uncertain. They are usually modeled starting from 
theoretical (e.g. [35, 36]) or observational (e.g. [37]) argu- 
ments. In particular, the same effect of increased power 



of angular correlations on large scales can be due to non- 
Gaussianity [33] or to a different model for the bias [37]. 
In order to make detailed predictions of the amount of 
non-Gaussianity that can be detected in forthcoming ra- 
dio surveys, it will be necessary to model the effects of 
variations of n(z) and b(z), but this is beyond the scope 
of this paper and is left to a future work. 



V. CONCLUSIONS 

We have presented an analysis of GR corrections to the 
angular power spectrum of radio sources that will be mea- 
sured with forthcoming radio continuum surveys such as 
LOFAR, EMU and WODAN, along with predictions for 
SKA-like surveys. These surveys will be well suited for 
probing GR corrections to the standard Newtonian anal- 
ysis of the spectra, since they will observe super-Hubblc 
scales not yet surveyed, being very deep and wide. 

We have computed for the first time the contributions 
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FIG. 6: The fractional difference of the galaxy count spectra between the cases with and without non-Gaussianity, Cf (/nl) — 
C<(/nl = 0)|/Cp R (/NL = 0). The solid and dashed curves illustrate the spectra with and without GR corrections respectively. 
Shaded regions show the cosmic variance. 



from rcdshift distortions and lcnsing convergence to the 
angular power spectrum for radio continuum surveys, 
thus generalizing previous results that incorporated only 
the overdensity contribution [4]. Then we have included 
the further GR corrections, in the form of velocity and 
potential terms. We have shown how all these contribu- 
tions will, individually and in combination, affect mea- 
surements of the angular power spectrum for the different 
surveys considered. The GR corrections to the standard 
Newtonian analysis are most significant on the largest 
scales, reaching 0(10%) for the SKA. However, precisely 
because they grow with scale, they are dominated by cos- 
mic variance for the near future SKA Pathfinder genera- 
tion surveys, while they could be in principle observable 
by a sensitive enough SKA-like survey. With complemen- 
tary information from other tracers (e.g. via the Euclid 
optical/ IR survey), cosmic variance can be overcome. 

A large-scale increase in power can also be due to pri- 
mordial non-Gaussianity. This means that a GR analysis 
is essential for a correct calculation of the non-Gaussian 
signal. We have computed the corrections arising from 
the local form of non-Gaussianity for the different sur- 
veys. Comparing the non-Gaussian effect to that of GR 
corrections without non-Gaussianity, we find that non- 
Gaussian corrections to the power spectrum will dom- 



inate over GR corrections for continuum surveys. The 
non-Gaussian signal rises above cosmic variance on large 
enough scales as follows: for WODAN when /nl ^ 80; 
for EMU when / NL > 50; for LOFAR when / NL > 20 
and for SKA when /nl ^ 5. 

Continuum radio surveys do not provide redshift infor- 
mation, so that GR corrections can be degenerate with a 
change in the distribution of matter, given by the prod- 
uct of radial source distributions with the bias. We ex- 
pect that an SKA HI galaxy survey will show stronger 
and more clearly defined GR corrections, because spec- 
troscopic information will break this degeneracy. 
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